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Abstract. We study Hamiltonian walks (HWs) on the family of three-dimensional 
modified Sierpinski gasket fractals, as a model for compact polymers in 
nonhomogeneous media in three dimensions. Each member of this fractal family is 
labeled with an integer b > 2. We apply an exact recursive method which allows for 
explicit enumeration of extremely long Hamiltonian walks of different types: closed 
and open, with end-points anywhere in the lattice, or with one or both ends fixed 
at the corner sites, as well as some Hamiltonian conformations consisting of two or 
three strands. Analyzing large sets of data obtained for b = 2, 3 and 4, we find 
that numbers Zm of Hamiltonian walks, on fractal lattice with N sites, for N ^> 1 
behave as Zn ~ w" fJ> N " ■ The leading term ui N is characterized by the value of the 
connectivity constant lu > 1, which depends on b, but not on the type of HW. In 
contrast to that, the stretched exponential term fi N depends on the type of HW 
through constant fi < 1, whereas exponent a is determined by b alone. For larger b 
values, using some general features of the applied recursive relations, without explicit 
enumeration of HWs, we argue that asymptotical behavior of Z^ should be the same, 
with a = ln3/ln[6(6+ 1)(6 + 2)/6], valid for all b > 2. This differs from the formulae 
obtained recently for Hamiltonian walks on other fractal lattices, as well as from 
the formula expected for homogeneous lattices. We discuss the possible origin and 
implications of such result. 

Keywords: solvable lattice models; polymers, copolymers, polyelectrolytes and 
biomolecular solutions; structures and conformations (theory) 

1. Introduction 

The fact that compact conformations of polymers are principal configurations of the 
native states of globular proteins makes them an important subset of all the physically 
accessible conformations. Since compact conformations occupy space as densely as 
possible, one of the simplest ways to model them is to use Hamiltonian walks (HWs) on 
a lattice, which are, by definition, self-avoiding walks (SAWs) that visit all the lattice 
sites exactly once [1]. In order to make this model more capable for capturing different 
features of various physical phenomena (such as protein melting [2] or protein folding [3]) 
local interactions can be introduced, but even in its simplest form, with no interactions 
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taken into account, the problem of enumeration and classification of Hamiltonian walks 
proved to be extremely difficult. 

Studies of Hamiltonian walks are primarily focused on finding the overall numbers 
Z N of open and closed HWs on lattice with N sites. It is expected that limiting value 
of In Zn/N exists when N — ► oo, and its particular value 

, ,. hiZ N 
muj = hm 

JV^oo N 

corresponds to the configurational entropy per monomer (site). This means that to 
the lowest approximation, Z^ behaves as u N and, therefore, the so-called connectivity 
constant uo > 1 can be interpreted as average number of steps available to the walker 
having already completed a large number of steps. The leading corrections to the 
exponential term are expected to have the power-law form N a (as in the case of ordinary 
self-avoiding walks), and the stretched exponential form fi N " , with fi < 1, so that Z^ 
should scale as 

Z N ~ u N n N "N a . (1.1) 

Expectations of such forms are based on the exact studies of HWs on the Manhattan [4] 
and some fractal lattices [5-7] , as well as on results obtained for closely related models 
of collapsed interacting self- avoiding walks on square [8-11] and cubic lattice [12]. In 
the case of collapsed SAWs appearance of the stretched exponential term, which is not 
present in the scaling form for non-interacting SAWs, was explained as a consequence 
of surface effects. Namely, a collapsed SAW forms a compact globule, with a sharp 
boundary separating it from the surrounding solvent, so that monomers on the boundary 
have smaller number of contacts with other monomers than those in the bulk of the 
globule, and therefore surface tension should arise. For homogeneous lattices one can 
assume that the boundary itself is homogeneous surface, and then straightforward 
arguments [13] lead to the conclusion that term /i^" 7 , with a = (d — l)/d (d being 
the dimensionality of the lattice), should appear. Exact studies [9], series analysis 
of data obtained via exact enumeration [8, 10], as well as Monte-Carlo simulations 
[11, 12] of collapsed SAWs on homogeneous lattices, certainly confirm existence of 
the stretched exponential term with the proposed formula for a. Although it is 
believed that HW model corresponds to the interacting SAW model at temperature 
T = 0, in spite of the continuous improvements of exact [14] and Monte Carlo [15] 
enumeration techniques, direct confirmation of scaling relation (11 .11) for HWs on non- 
oriented homogeneous lattices has not been achieved yet. Apart from Manhattan lattice, 
scaling forms for Hamiltonian walks were obtained only for some fractal lattices, where 
stretched exponential corrections for both open and closed HWs were found only for 
ra-simplex fractals with even n [7] . Even in that cases, its presence cannot be explained 
using a simple generalization of the argument used for collapsed self-avoiding walks 
on homogeneous lattices. However, results of the studies on fractals suggest that 
stretched exponential term can be expected on lattices where larger number of entangled 
conformations is possible. In order to get a deeper insight into this issue, here we apply 
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an exact recursive method for enumeration and classification of HWs on modified three- 
dimensional Sierpinski gasket family of fractals. In contrast to most of the previously 
in this context studied fractals, which were embedded either in two-dimensional spaces 
or spaces with dimensionality higher than 3, each member of the family studied in the 
present paper is embedded in three-dimensional space. Therefore, HWs on these fractals 
can be understood as a toy model for compact polymers critical behaviour in realistic 
nonhomogeneous 3d media, which, to the best of our knowledge, has not been studied 
so far. 

The paper is organized as follows. Modified 3d Sierpinski gasket fractals, as 
well as the method itself, are described in the next section. Explicit forms of the 
recursion relations, obtained for the types of Hamiltonian walks which are needed for 
generation and enumeration of all closed conformations are presented for the first three 
members of the fractal family. Numerically analyzing these relations we find scaling 



form a; . Using some general features of these relations (obtained in Appendix A 



and Appendix B ), we argue that such form should be correct for the whole fractal 
family, and derive the closed-form formula for the exponent a. In section 3 we extend 
the method to open HWs, explicitly apply it again on the first three fractals, and then 
generalize it. It turns out that the number of open HWs scales as uo N ^,q" , with the same 
values of the connectivity constant u and the exponent a as for the closed HWs, but 
with fio 7^ f'C- General features of the recursion relations for open HWs are derived in 



Appendix C| and |Appendix D| All the obtained results are summarized and discussed 



in the section 4, and some auxiliary considerations, needed for the discussion, are given 
in Appendix E. 



2. Closed Hamiltonian walks on modified three-dimensional Sierpinski 
gaskets 

Three-dimensional modified Sierpinski gasket (3d MSG) fractal is constructed 
recursively, starting with a unit tetrahedron. The first step of the construction, so 
called generator Gi(b) of order / = 1, is obtained by joining 

N G = b(b + l)(b + 2)/Q (2.1) 

unit tetrahedrons into a b times larger tetrahedral structure (see figure [1]) in such a way 
that vertices of neighboring unit tetrahedrons are connected via infinitesimal junctions. 
Enlarging the generator b times, and substituting the smallest tetrahedrons with G\(b), 
and then repeating this procedure / times, one obtains Gi + \(b) - generator of order 1+ 1, 
which contains Ni + i = ANq' 1 sites. The complete 3d MSG fractal with parameter b is 
obtained when I — > oo, and its fractal dimension df is equal to 

d / = ln[6(6 + l)(6 + 2)/6]/ln6. (2.2) 

One should note here that ordinary 3d Sierpinski gasket (SG) fractal with parameter 
b is constructed in a similar way, with the only difference being that neighboring unit 
tetrahedrons in 3d SG lattice are not moved apart, as is the case for the 3d MSG fractal. 
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Figure 1. The first two steps of the construction of the 3d MSG fractal with b = 4. 
Vertices of neighboring unit tetrahedrons are connected via infinitesimal junctions. 



This small change does not alter the basic properties of the lattice, such as its fractal 
dimension, but it simplifies the complicated scheme needed for enumeration of closed 
Hamiltonian walks on 3d SG fractals [6] and allows for enumeration of extremely long 
open HWs on 3d MSG fractals with 6 = 2,3, and 4. In addition, general conclusions 
about scaling forms of the numbers of HWs on 3d MSG fractal with arbitrary b > 2 can 
be derived, as will be explained in what follows. 

By definition, each Hamiltonian walk on Gi(b) structure visits all its Ni sites exactly 
once. In figure[2]an example of closed Hamiltonian walk on the generator G\ (3) is shown. 
One can notice that this HW can be decomposed into Nq = 10 parts (within the ten 
unit tetrahedrons), which are either one- or two-stranded. This observation can easily 
be generalized to generators of higher order and any b: any closed HW on Gi + i(b) can be 
decomposed into Nq parts within the same number of generators Gi(b), from which the 
generator Gi + \(b) is made of. These parts traverse each Gi(b) either once or two times, 
and we shall call the corresponding configurations B- or E-type HW steps, respectively, 
whereas numbers of HWs of these types within a G\(b) will be denoted by Bi and Ei, 
respectively. Using these numbers, the overall number of closed HWs on Gi + %(b) can be 
expressed as 

Z? +l = Y.n k B^- k E^ (2.3) 

fc=0 

where n k are the numbers of closed HWs configurations within Gi+i(b), with k steps 
of E—type, and (Nq — k) B-steps. For instance, the path presented in figure [2] is one 
of Zf possible closed HWs within the G\(3), contributing to the term u 2 BqEq in the 
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Figure 2. Example of closed Hamiltonian walk (red line) on the generator of the 3d b = 
3 MSG fractal. The generator consists of 40 vertices, which are numbered consecutively 
from the first (arbitrarily chosen) vertex visited by the HW, to the last one (which is 
connected with the first one). One can notice that HW conformation within the unit 
tetrahedron can be either one-stranded (B-type), like the one connecting the sites 26- 
27-28-29, or two-stranded (£7- type), like the conformation traversing the tetrahedron 
with vertices 4-5-12-13. 



corresponding equation ( 12. 31) . The upper limit kc in the sum in ( 12. 31) is equal to for 
6 = 2, whereas for b > 2 it can be shown (see Appendix A ) that 

k c = ^(6+l)(6 + 2)-8. (2.4) 

Due to the self-similarity of the lattices under study, the numbers n k do not depend on 

I, and, also, numbers B\ and E\ fulfil recursion relations of the following form: 

N G N G 
B l+l = £ m k B k E^~ k , E l+1 = £ p k B k E^ k , (2.5) 

k=kg k=kg 

where coefficients m k and p k depend only on k and b, and are positive integers for all 
b > 2, including zero for b = 2. For instance, for b = 2 these relations ar^j] 

B l+1 = 2Bf + ABfEi + 6B?E? , E l+1 = Bf + ABfE l + 22E? , (2.6) 

and we were able to find the explicit form of relations (12.51) for b = 3 and 4 also, by 
direct computer enumeration of possible HW conformations within the MSG generator. 
The corresponding coefficients m k and p k are presented in table [H while for larger b 
values they could not be reached within the reasonable time with the computer facilities 

| Note that for 6 = 2 recursion relations ()2.5p are the same as for the 4-simplex fractal lattice [5] . 
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Table 1. Coefficients appearing in recursion relations (|2.5|) . found by direct 
computer enumeration of the corresponding HW conformations on 3d MSG fractals 
with 6 = 2,3,4. 
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available to us. However, one can show (see Appendix A ) that for general b > 2 lower 
limits of the sums in relations (12.51) are equal to 

k B = N G - k c - 3 = -{b + 1)(6 + 2) (6 - 3) + 5 , k E = k B - 3 , (2.7) 

o 

which is a result essential for establishing some general conclusions, as will be explained 
in the following paragraphes. 

Once the explicit form of recursion relations ( 12.31) and (12.51) is established, starting 
with the initial values Bq = 2, and Eq — 1, corresponding to the unit tetrahedron, one 
can calculate the numbers B h E h and subsequently the overall number Z^ +1 of closed 
HWs, in principle for any /. However, since these numbers quickly become extremely 
large, it is useful to introduce new variable x\ = B\jE\. Then, from (I2.5P it follows that 
recursion relation for the variable x„ for b > fi has the form 

fcc+3 

E m kB+k x k 

x l+1 = x^g^) = xf (/(i,) , g{x) = ^ . (2.8) 

E Pk E +kX k 

k=0 

Numerical analysis of this recursion relation for b = 3 and 6 = 4 cases, shows that 
starting with the initial value xq = 2, x\ quickly tends to 0, as / grows. Assuming 
that this is correct for general b > 2, for large values of / the above equation can be 
approximated as 

x l+l ^^xf, (2.9) 
Vk E 

§ The corresponding analysis of the 6 = 2 case is the same as for 4-simplex fractal lattice, which is 
given in [5]. 
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Table 2. Values of relevant constants appearing in the scaling forms (|2.17|) and (|3. 13|) 
of the overall numbers of HWs on 3d MSG fractals, with 6 = 2,3, 4, together with the 
corresponding values of the fractal dimension d / . 
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from which follows 

xi ~ A 3 ' , (2.10) 

where A is some constant, whose value depends on b, but it is always less than 1. For 
b = 2 the corresponding relation is x\ ~ A 2 ' [5], and particular values of A for b = 2,3 
and 4 cases are given in table El 

Recursion relation for numbers E\ of two-stranded HWs (12. 5p . with the variable 
x\ = Bi/Ei obtains the form 

fcc+6 

E l+1 = Efoaf'fixi) , f(x) = Pk + k E x k , (2.11) 

k=0 

from which one gets 

lnK + i InEi , lnxz In f (xi) 

ttt = r + k E tttH -Att-- (2.12) 

Numerically iterating this recursion relation, together with ( 12. 8ft . one finds that 

lim — —r- = In uj , (2.13) 

i^oo 4N l G K ' 

where uj (see table [2]) is constant larger than 1. On the other hand, the overall number 
Z^ +l ( 12.31) of closed HWs on the generator of order I + 1 can be expressed as 

Zf +1 = Ef G x^- kc h{ Xl ) , h(x) = £ n kc . k x k , (2.14) 

fc=0 

so that 

\n Zf +1 _ In Ei N G -k c \nxi \nh(x t ) 

AN£ l ~ AN l G AN G N l G AN£ l ' 1 ' ' 

From the asymptotical behavior (12. 10)) of the number xi, and from the fact that h(x) 

tends to constant value n kc > when x — > 0, it then follows that 

In Zf \nEt , 
lim = lim j- = In uj , 

i^oo Ni i^oo AN l G 

where Ni = AN l G is overall number of vertices within the generator of order I. 

To find the leading-order correction to the asymptotic behavior of Zf we first 
introduce variable 

\nE t lnE t 

Vl= AN G - = ^' (2 - 16) 
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which, as follows from (12.121) . satisfies the relation 
i-i 1-1 i 

Vi = XX^'+i -Vk)+yo = J2 — [kElnx k + \nf(x k )] . 

k=0 k=0 ^k+l 

Then, using (12.131) . one obtains 

OO 

yi = \nuj - T7 — [k E \nx k + \nf(x k )} , 

k=l ^k+l 

from which, taking into account the large k behavior (12.101) of x k , it follows that 

i 1 o/ , , const , 
y t « In u — 3' In A , I > 1 . 

Substituting this relation into (12.151) one derives 

\nZf « Nilnuj + Nflnfic, i.e. Zf^uj^^J, (2.17) 

with 



and 



m 3 In 3 , . 

° = LTiVc = ln[6(6 + 1)(6 + 2)/6] ' ( ' ' 

,c = X\ A=^±^A^. (2.19) 

At the end of this section we want to emphasize that the crucial step in preceding 
derivation, which led to the scaling form ( 12.171) . with the value of exponent a given by 
formula (I2.18p . was the assumption that xi tends to for general b (explicitly confirmed 
only up to b = 4). Its direct consequence is approximate difference equation (I2.9p . which 
is obtained from the exact relation (12.81) . On the other hand, the key ingredient of that 
relation is the fact that ks — k E = 3, due to which xi behaves as A 3 ', for I 3> 1. In that 
sense, exact expressions for ks and k E ( 12.7[) . obtained in Appendix A for general b > 2, 
are essential for establishing the scaling form (12.171) . whereas the particular values of the 
coefficients m k and p k of the recursion relations (12.51) do not affect neither its general 
form, nor the value of a. However, whether x\ tends to or not certainly depends 
on the values of m k and p k . Analyzing data given in table [1] for b = 3 and 6 = 4, 
one can observe that p kE + k > m kE+k+ i, for 2 < k < k c + 5. It can be shown (see 



Appendix B) that such inequality for general b > 2 is sufficient to prove that numbers 



xi, which satisfy difference equation ( 12.81) with the initial condition xq = 2, tend to 
when / — > oo. Unfortunately, the inequality itself we were not able to prove, but, since 
it seems plausible, we think that scaling relation ( 12.171) can be accepted as valid for 
general b. 



3. Open Hamiltonian walks on modified three-dimensional Sierpinski 
gaskets 

Any open HW on generator Gi+\{b) of order / + 1 can be decomposed into Nq parts 
within its iV^ constitutive generators Gi(b) of order /. The parts which contain ending 
points (see figure [3]) can be of four different types: 
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• A-type, which consists of one HW, with one end at one vertex of Gi(b) and the 
other end at any other site of Gi(b), including its vertices; 

• C-type, composed of two non-intersecting SAWs, one with both ends at the vertices 
of the Gi(b), and the other with one end at the third vertex of the Gi(b) and the 
second anywhere within it - these two SAWs together visit all the sites of the Gi(b); 

• D-type, consisting of two non-intersecting SAWs, each of them starting at different 
vertex of Gi(b) and ending anywhere within it, in such a way that all sites of Gi(b) 
are visited; 

• H-type, comprised of three non-intersecting SAWs, each of them starting at 
different vertex of Gi(b), one ending at the fourth vertex, and the remaining two 
anywhere within Giib) - again, all sites of Gi(b) should be visited by these three 
SAWs. 

A- and C-type configurations have one dangling end, therefore we shall call them one- 
leg configurations (steps), whereas D and if -configurations, with two dangling ends, 
will be called two-leg steps. If both end-points of the complete HW lie in the same 
Giip), that Gi{b) contains a two-leg step (of D- or if-type), whereas the remaining 
Nq — 1 generators Gi(b) contain either B— or Estep. Otherwise, when end-points are 
in two different Gi(b) generators, the complete open HW has two one-leg steps (of A— 
or C-type), and the remaining Nq — 2 parts of the walk are of either B- or J5— type. 
Therefore, one concludes that overall number Z° +1 of open HWs within Gi + \(b) is equal 
to 

Z° +1 = AfF AA + AiCiF AC + CfFcc + DiS D + HiS H , (3.1) 

where F AA , F A c, Fee, Sd and Sh are polynomials in Bi and Ei, of power Nq — 2 
or Nq — 1, and A h Ci, D\ and Hi are overall numbers of corresponding type HWs 
within Giib) generator. Numbers of one-leg conformations fulfill recursion relations of 
the following form: 

Ai+i = Rn(Bi, Ei)Ai + Ri2{Bi, E{)Ci , ^ 

Ci + \ = R 2 i(Bi, Ei)Ai + R 2 2(Bi, Ei)Ci , 

with 

R ti (B,E) = irr%B N °- l - k E*, (3.3) 

k=0 

while numbers of two-leg conformations obey relations of the form 
A+i = d D Di + d H Hi + d AA Af + d AC AiCi + d cc C? , 
H i+1 = h D Di + h H Hi + h AA A] + h AC AiCi + h cc Cf , °'" > 

with d and h being some polynomials in Bi and E\. For instance, for 6 = 2 relations 

(ED) and (E2D are: 

Zf> +1 = 12Bf(Af + 2AiCi + 3Cf + BiDi) , 

A t+1 = (6B? + QBfEi)Ai + (12 B? + 18flf , 

C i+1 = (Bf + 3BfEi)Ai + (3Bf + UBjEi + lQBiEf + 16Ef)Q . 
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Figure 3. Three examples of open HW (oriented red lines, where the only meaning of 
arrows is to serve as guides to the eye) on the generator of the 3d b — 3 MSG fractal. 
Points at which HW begins or terminates are marked as full red circles, for the sake 
of better recognition. If both ends of the walk are in the same tetrahedron the two- 
leg HW conformation within that tetrahedron is either of D- or H-type. Otherwise, 
if end-points of the walk belong to different tetrahedrons, the corresponding one-leg 
conformations are of A- or C-type. 



They coincide with the corresponding relations for 4-simplex fractal lattice, as well as 
does the complete further analysis, which can be seen in [7]. In the remaining part of 
this section we present the general analysis for b > 2 3d MSG fractals. 



For general b > 2 it can be shown (see |Appendix C\j that polynomials F and S, 
appearing in (13.11) . have the form 

k c +i 

F XY (B, E) = z^ Y B N °- k ' 2 E k , XY = A A, AC, CC 

S D (B, = £ zj? B Na ~ k ~ 1 E k , S H {B, E)= £ z^B^'^E* . 

k=0 k=0 

Starting with the initial values for the numbers Bi, Ei, A[, Ci, Di, and Hi, and using the 
recursive relations (12. 5p . (13. 2p . and (13.41) . one can calculate these numbers in principle 
for any I, and substituting them into (13. If) . eventually find the overall number Z^ +l of 
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open HWs. Since all these numbers increase rapidly with I, it is useful to introduce 
variables 

A G x Di H t 

ui = —, Vi = -=p , Wi = -=- , qi = — , (3.6) 
tbi tii tii tii 

in addition to already defined xi = Bi/Ei. With these variables ZP 1 can be rewritten 

as 

'k +l 

Zf \ = E^ G x l NG ~ kc ~ 



'7 1 X k 



k=0 



kc kc—1 

+ nx 2 i E z k C -k x i + n x i E z k C -k-i 

k=0 k=0 



(3.7) 

and explicitly calculated using corresponding recursion relations for ui, vi, wi, and qi, 
simultaneously with ( 12. 8ft and (12. lip . In Appendix C it is shown that the upper limits 
of the sums in (13.31) are equal to 

hi = hi = kc + 2 , hi = k 22 = k c + 5 , (3.8) 

so that by dividing ( 13.21) with recursion relation for Ei, given in ( 12.51) . for ui and vi one 
obtains recursion relations 



(3.9) 



/ ui + i \ _ I mn | «i 

\ v l+1 J \ m 2 i m 22 J \ vi J 

where rriij are functions of Xi of the following form: 

x 3 k c +2 j k c +5 

Tn>ii(x) = 77 \ E r fcc+2-fc x > m 2i(x) = ^ r fe ^ +5 _ fe ;r , z = 1,2, 
A X J fc=o A X J fe=o 

and /(x) is defined in (12.111) . Since £/ — > and /(^) — > const ^ when Z — ► oo, it 

follows that m\i{xi) — > 0, m 2 i(xi) — » const 7^ 0, implying that tends to 0, and to 

some constant value. This is indeed correct for 6 = 3 and 6 = 4 cases, for which we 

managed to find the complete set of coefficient r k J (see table [3]). Starting with the initial 

values A = 6, Co = 2, i.e. uo = 6 and vq = 2, and numerically iterating relations ( 13. 9ft , 

for both 6 = 3 and 6 = 4, already after five iterations one obtains limiting values of vf. 

441.32 ... for 6 = 3, and 3538.91 ... for 6 = 4. 



In Appendix D it is shown that recursion relations for the variables wi and qi can 
be put into the following matrix form 

( Wl+1 \ = ( Pn Pl2 \ ( Wl \ + ( tnU * + tl2UlVl + tl3V ? \ (3 10) 
\ Ql+i J V P21 P22 ) \ Ql ) V *2iU? + t^uivi + t 23 vf J 

where pij and tij are functions of xi of the form 

x 2 fcc+3 2x 3 kc+2 

Pn(x) = — E ( fc + k B )m k+kB x k , pi 2 (a;) = — — (fee + 3 - k)m k+kB x k , 
k =o J \ x ) k =o 

P2i(a?) = 9 77 \ E ( fc + k E )Pk+k E x k , P22W = 77-^ E (*o + 6 - k)p k+kE x k , 
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Table 3. Coefficients of the polynomials Rij (|3.3|) . appearing in recursion relations 
(|3.2j) . found by direct computer enumeration of the corresponding HW conformations 
on 3d MSG fractals with 6 = 3 and 4. 

6=3 6=4 



k 


' k 


' k 


r 21 


' k 


' k 


' k 


r 21 
' k 


r 22 
' k 





378 


1050 


202 


620 


726000 


2078460 


454340 


1343670 


1 


1722 


5256 


798 


2854 


6785652 


20162052 


4124238 


12825924 


2 


2748 


8964 


1962 


8696 


30159948 


92803776 


18522682 


61201186 


3 


3588 


15144 


2912 


15172 


85107696 


270289488 


56914854 


200289608 


4 


2976 


9576 


3126 


24064 


171919920 


559983864 


135608360 


503818972 


5 






4632 


23684 


265804824 


875876784 


260879696 


1011455276 


G 






2420 


27088 


319608912 


1041132024 


412200376 


1653387800 


7 






1304 


4308 


290292720 


895898352 


548029844 


2243585880 


8 










185693208 


493595664 


612939852 


2479455304 


9 










51662736 


71546160 


552106184 


1451395108 


10 














407928692 


1451395108 


11 














226652724 


546796688 


12 














68954256 


26465392 



tidx) 



fc c +4 

E ■ 



XY 
'fc c +4-fc 



X 



t2l(x) 



k c +7 

E « 



k=0 



-.AC 



E °kc+6-k X 



x 3 f(x) ^ 

1 fec+5 
fc=0 



AA 

kc+7-k 



X 



cc 

'k c +5- 



,X 



(3-112 



z 2 /(aO fc=0 

with and o^ y being some positive constant integers. Using these relations, one can 
calculate W[ and qi, and by putting them, together with xi, ui, vi and E h into (13. 7p . 
finally can evaluate Z®, for any /. 

In order to find asymptotic behavior of Z®, we first notice that relations (I3.10p can 
be combined with the relation ( 12. 8ft for xi, so that one can establish recursion relations 
for Wixf and qixf, terms through which variables Wi and qi appear in (13.71) . Using the 
facts that x\ — > 0, f{x{) — > pk E , Ui — > 0, and V\ — > const, when I — > oo, it can be shown 
that 





(3.12) 



where a,;,- and 613 are some constants (see 



Appendix D). This relation implies that wixf 
00, and consequently, comparing the terms containing Wixf 



and qixf tend to zero as I 
and qixf with those with ui and in (13.71) . one can conclude that 



7° 



pN G Na-kc- 
x l 



fc=0 



k c +l-k a l + ^ 



This means that for / ^> 1 number behaves as 



7 o 



3 



N G N G -k c - 
x l 
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so that, using (I2.14p . one obtains 
Z ?+i -3 \-3 !+1 

Consequently, taking into account (j2.17p . it follows that overall number Z° of open 
HWs scales as 

Zp~u" l ltf, »o = \ B , B = ^ ± ^^ ( 3 - 13 ) 

where, as in (I2.18p . a = ln3/lniVG, and particular values of u, fio an d ^ are given in 
table [2J 



4. Summary and Discussion 



In this paper we have analyzed asymptotic behavior of the numbers of open and closed 
Hamiltonian walks on three-dimensional modified Sierpinski gasket family of fractals. 
Numbers of extremely long HWs on these lattices can be generated by applying an exact 
recursive enumeration scheme, based on the fact that any HW on the (I + l)th step of 
the construction of the 3d MSG fractal, Gi + i(b), consists of No = b(b + l)(b + 2)/6 
HWs within its Nq constitutive Gi(b) structures, which can be of one of six possible 
types: A, B,C, D, E, and H (figure [2] and figure [3]). Numbers of these HWs, A h B h 
Ci, Di, Ei, and Hi, fulfil closed set of recursive relations, (12.5j) . (13. 2ft and (13.4ft . which, 
due to the self-similarity of the lattices under study, do not depend on I, and therefore 
can be obtained by explicit enumeration and classification of HW conformations on the 
first step of the fractal construction, G\(b). As b grows, number of HW conformations 
rapidly increases already on G\(b), so that we have managed to find explicit form of the 
recursive relations only up to b = 4. However, we have shown that for any b > 2 these 
recursive relations have some features that enable general analysis, leading to conclusion 
that overall numbers Zf and Zf of closed and open HWs, respectively, on any 3d MSG 
fractal, scale with the number Ni 



AN G of the lattice sites as 



A" 



Ng+ k C . 



with cr given by formula ( 12. 18ft : 

In 3 



In 3 1 



k C+ 3 A- 

\ N G~ 3 



Nf 



a 



ln|6(6 + l)(6 + 2) lnbd f ' 
Constants u and A can be obtained numerically, using relations 



lim 

l— >oo 



InK 



In A 



lim 

Ni i->oc 3 l 

and their particular values for b = 2, 3, and 4 are given in table [2j Number 
kc = {b + 1)(6+ 2)/2 — 8 represents the maximal number of Gi(b) generators within the 
Gi+i(b) they belong to, which are traversed by two-stranded E-type HW conformation 



HBi/Ei 



(E-step) within any closed HW on Gi + \{b) (see equation (12. 3p and | Appendix Afl . As 
one can see, HW conformations that traverse Gi(b), and recursion relations (12.51) for 



A model of compact polymers on a family of three-dimensional fractal lattices 14 

the corresponding numbers B t and E t , are sufficient for obtaining connectivity constant 
uj that governs the leading exponential term in the scaling forms for Zf and , as 
well as constant A which appears in their correction terms. Whereas the leading term 
u Nl is the same for both closed and open HWs, correction term in both cases has the 
same stretched exponential form /i N ^ , but with different values for /i and \iq > /ic, 
indicating that number of open HWs is larger than the number of closed HWs. One 
should notice here that the fact that scaling form obtained for open HWs is determined 
by the behavior of the numbers Bi and E\ only, means that contribution of one- and 
two-leg HW conformations (A-, C-, D and if-type walks), i.e. HWs with their ends in 
the interior of Gi(b) fractal structures, is not significant. Furthermore, in sections 2 and 
3 it was elaborated that the main mathematical reason for obtaining scaling forms with 
the stretched exponential correction term is the fact that Mmi-t^Bi/Ei = 0, meaning 
that entangled conformations (with large number of two-stranded parts) dominate over 
those in which HWs rarely return to already visited fractal generators. 

Results of the study of HWs on 3d MSG fractals, presented in this paper, should 
be compared with the results recently obtained for other fractal lattices [7]. For n- 
simplex fractal lattices with even values of n the same asymptotic behavior of the 
HW numbers was found, uo N [i N " , but with a = l/df, which differs from the formula 
derived here for a. The existence of the term was also related with the facts that: 
(1) entangled HWs prevail, and (2) HWs with "interior" terminating points do not 
contribute to the asymptotic behavior of HW numbers, as in the 3d MSG fractals case. 
For odd n, stretched exponential terms in the scaling forms for HWs on n-simplex 
fractals do not exist, HWs with "interior" terminating points do affect the scaling 
forms for the numbers of open HWs, and entangled conformations do not dominate 
in this case. For Given-Mandelbrot and 2d modified Sierpinski gasket fractals (which 
are two-dimensional generalizations of the Sierpinski gasket and 3-simplex fractal, 
respectively) stretched exponential terms are also absent, whereas again, one- and two- 
leg conformations are necessary for obtaining the scaling forms for open HWs, and, 
due to the specific topology of these lattices, only one-stranded HW conformations are 
possible. Therefore, results of the present study confirm the assumption that existence of 
the stretched exponential term is related with the issue of HW entanglement. However, 
the discrepancy between the formulas obtained for the exponent a for different fractal 
families, and, more generally, the question which properties of the underlying lattice 
and in what way determine a, deserves further consideration. To this end we recall that 
stretched exponential terms were obtained for HWs on Manhattan [4], as well as for low- 
temperature SAWs on square [8-11] and cubic [12] lattice. In this context, one should 
also mention the spiral SAWs on 2d regular lattices [16] and lattice animals on some 
hierarchical lattices. Whereas, to the best of our knowledge, physical interpretation of 
the stretched exponential term for the 2d spiral SAW models has never been proposed, 
its presence in the scaling forms obtained for lattice animals on hierarchical lattices 
was directly connected with the existence of the sets of sites with different coordination 
numbers [17]. Utilizing similar idea, in the case of low-temperature (collapsed) SAWs it 
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was explicitly explained in [13] why term of the form fJ, N<r , with a = (d — l)/d for regular 
lattices, should arise. However, as was elaborated in [6,7], direct application of such 
approach on HWs on fractal lattices does not give satisfactory result. Yet, the scaling 
forms obtained in this paper can be expound in the spirit of the physical reasoning given 
in [13], as will be explained in the following paragraph. 

We shall first focus on closed HWs on Gi(b) with the maximal number of two- 
stranded E conformations within the unit tetrahedrons. Such HWs represent maximally 
entangled closed (MEC) compact conformations and they accomplish the maximal 
possible connectedness between the generators Giip) of all orders /. Now, from all the 
Ni sites visited by such maximally entangled HW, observe those which belong to one- 
stranded conformations (5-steps) within the unit tetrahedrons, and which are directly 
connected only with the sites belonging to the same unit tetrahedron. For instance, 
vertices 9 and 10, or 15 and 16 in figure [2] represent examples of such sites. Since 
these sites are not directly connected with the other tetrahedrons of the lattice (via 
HW), they are the maximally isolated sites within the maximally entangled closed HW 
conformation. It is shown in Appendix E that number Nf 1 of such sites is equal to 

1 2 A^ G -3 1 N G -Z 1 ' V ; 

and one can easily see that, using this expression, scaling relation for the numbers od 
closed HWs can be transformed into the following form 

/—N CI N G- k C~ 6 

Zf~uj' Nl V\ l , with J = lj/X 4 ^g- 3 > . (4.2) 



In a similar way (see Appendix E ), one can show that scaling relation for open HWs 
can be expressed as 

Z° ~u' Nl V\ N?1 , (4.3) 

where 

' 2 Nq — 3 Nq — 3 1 K J 

is the number of maximally isolated sites on Gi(b) visited by maximally entangled open 
HW (which is an open HW with the maximal number of .E-steps on unit tetrahedrons 
and with both loose ends being of type C, on all levels I). Comparing forms ( 14.21) and 
( 14. 3p with the scaling relations expected for HWs on homogeneous lattices, one can say 

/—N 01 j—N 01 

that in this case terms v A 1 and v A ' play the role of the stretched exponential 
correction term fi N ° in the case of homogeneous lattices. Indeed, for homogeneous 
lattices N a , with a = (d — l)/d, is proportional to the number of sites on the lattice 
boundary, which have smaller number of neighbors than the bulk sites, and in that sense 
they are similar to the maximally isolated sites within the maximally entangled HWs 
on 3d MSG lattices. Therefore, one might generally expect that scaling relation for the 
number of HWs on nonhomogeneous lattices instead of stretched exponential correction 
term has a term fi Nl , where Nj is the number of conveniently defined "maximally isolated 
sites" . Relation between Nj and the overall number of lattice sites A^ depends on the 
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particular topology of the lattice under study, which is the reason why different scaling 
forms, as functions of N, are obtained. 

Conclusion of the previous paragraph can be supported by performing a similar 
analysis in the case of HWs on previously studied fractal families [7] . For instance, for 
ra-simplex fractals maximally isolated sites can be defined in a similar way as for 3d MSG 
fractals, i. e. these are the sites which are directly connected by maximally entangled HW 
only with the sites within the same fractal generator they belong to. Here, maximally 
entangled HWs are those that contain maximal number of maximally-stranded HW 
parts. Then, it can be shown [18] that Nf 1 = 2JVjf and N? 1 = (f - l) Nf, where 
a = In 2/ Inn. Since overall numbers of closed and open HWs on n-simplex fractal with 
even n asymptotically behave as Zf ~ u Nl \ B 1 (see (5.8) in [7]) and Zf ~ u Nl \ B 2 1 
(formula (5.14) in [7]), respectively, they immediately fit into the scaling relation uo' N fj, 1 , 
where uj' = u and fi = Xb for both types of HWs. For odd values of n, however, it 
can be shown [18] that number Nj of maximally isolated sites, for both open and closed 
HWs, is proportional to iVj, so that any term of the form fj, Nl can only contribute to the 
leading term uj n . This certainly is in accord with the scaling forms found in [7] for odd 
n, which do not have stretched exponential corrections. 

To conclude, we may say that presented exact study of Hamiltonian walks on 
modified three-dimensional Sierpinski gasket fractals supports the idea that entangled 
conformations are of the most physical importance for the behavior of compact polymers 
in inhomogeneous media. The introduced concept of "maximally isolated sites" 
proved to be useful in physical interpretation of different scaling forms obtained for 
Hamiltonian walks on all so far studied fractals. Whether such concept can be applied 
on homogeneous and other nonhomogeneous lattices will be the matter of our further 
investigations. 
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Appendix A. Maximal number of E'-steps within HW configurations 

In this Appendix we prove that numbers kc, ks and k,E are given by formulae ( 12.41 ) 
and (12. 7p . The number kc is the maximum number of -E-steps within the closed HW 
configuration (see (12.31) ). whereas the numbers k B and k B are related to the maximum 
number of E-steps within the HW configurations of B and E-type (see (12.51) ). i.e. the 
conformations that traverse the generator of order I + 1 once or two times, respectively. 
By ,,i?-step" here we imply any HW conformation consisting of two mutually avoiding 
strands, both traversing the generator of order I (see figure[2]). In a similar way, ,,-B-step" 
is any HW conformation that consists of one strand traversing Gi(b). 

According to our definition of the modified Sierpinski gasket fractals, generators 
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Figure Al. On the left-hand side of this picture generator G;+i(&) of order I + 1, for 
b = 4 3d MSG, is presented, as it is seen from above. Gray-shaded small tetrahedrons 
represent generators G;(4) of order /, lower layers being darker. Vertices lying in the 
same horizontal plane are indicated with the same symbol, being also the same for two 
adjacent planes containing vertices that belong to different generators. On the right- 
hand side of the picture horizontal layers of tetrahedrons are split and slightly magnified 
in order that all vertices, as well as some junctions can be seen. In particular, twofold 
junctions in the lowest layer of vertices are indicated by yellow lines, whereas junctions 
joining three vertices are given in green color. There is only one four-fold junction 
- the corresponding four connected vertices are encircled, and in the magnified circle 
these vertices are presented as seen in three dimensions, with the junction indicated 
in orange color. One can check that overall number of junctions is Nj — 31 (three 
twofold junctions for each of six G;+i(4) edges, three three-fold junctions for each of 
four G;+i(4) faces, and one four- fold interior junction), which certainly is in accord 
with formula (jA.lj) . 



Gi(b) within the generator Gi + \(b) are connected via infinitesimal junctions. Each of 
these junctions connects two, three or four neighboring Gi(b) generators. ,, Twofold" 
junction connects vertices of neighboring Gi(b), both lying in the same edge of Gi + i(b) 
(for instance, vertices 3 and 4 in figure [2] are joined by such junction). Vertices of 
neighboring Giib) that lie inside the faces of Gi + i(b) (such as sites 5,6 and 19 in figure [2]) 
are connected by three-folded junction, whereas four-folded junction occur in the interior 
of Gi + i(b), and can exist only for b > 4. In figure [ATI we explicitly indicate some of the 
junctions connecting the vertices of the generators Gz(4), within the generator Gi + \(4). 
It is not difficult to show that the number of junctions which connect Nq generators 
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Giib) into the generator Gi + \(b) is equal to 

jV J =i(6 + l)(& + 2)(& + 3)-4. (A.l) 
o 

From the restraint that each site of the lattice has to be visited exactly once, follows 
that each junction (i.e. its middle point) can be traversed at most once. Now, suppose 
that closed HW conformation within Gi + i (b) consists of a £>-steps and (3 .E-steps. Each 
Giib) is traversed by HW, which implies that 

a + f3 = N G = h(b+l)(b + 2). (A.2) 
o 

Each £?-step uses two junctions and each E-step uses four junctions. However, since 
each junction connects two Gi{b) generators, the overall number of visited junctions is 
equal to a + 2/5. In addition, steps through corner Gi(b) generators for the closed HW 
have to be of the B-type (see figure [2]), which means that one of the three junctions for 
each of four corner generators is certainly not used. This implies that inequality 

a + 2(3 < Nj - 4 = \{b + 1)(6 + 2) (6 + 3) - 8 (A.3) 
o 

is satisfied. Then, from (IA.2I) and flA.3j) directly follows that 

(3<^(b+l)(b + 2)-8, (A.4) 

i.e. the maximum number of .E-steps within the closed HW conformation is indeed 
equal to 

Anax = k c = \{b + 1)(6 + 2) - 8 . (A.5) 

Next, we consider a £>-type HW configuration on generator, with a B- 

steps, and j3 E-steps. Steps through two corner Gi(b) generators, at which HW starts 
or terminates, can be of either B— or E— type. However, both remaining two corner 
Giib) generators must be traversed by 5-step (see figure IA"2| . meaning that two of Nj 
junctions are never used. On the other hand, each I?-step uses two junctions, unless it 
is the first or the last step of the HW configuration, in which case it utilizes only one 
junction (since the corner vertex does not belong to any junction). Similarly, E-step 
which traverses interior Giib) generators uses four junctions, whereas possible i?-step 
through either the first or the last traversed Gi(b) uses only three junctions. Suppose 
that both the first and the last visited corner generators contain E-step. Then, these two 
steps together use six junctions, and the remaining (J3— 2) ,, interior" E-steps use 4(/3 — 2), 
whereas a .B-steps use 2a junctions. Since every junction can be traversed only once, 
and two consecutive steps share junctions, it follows that altogether [6 + 4(/5 — 2) + 2a]/2 
junctions are visited by such HW conformation. This implies the following inequality 

^[6 + 4(/3-2) + 2a]<iV J -2, (A.6) 

which, together with the relation a + (3 = Nq, gives 

(3 < Nj - N G - 1 = i(6 + 1)(6 + 2) - 5 = /3 max = N G - k B . (A.7) 
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Figure A2. Generator G/+i(5), viewed from above, with the upper three layers of 
constitutive generators Gz(5) (gray-shaded tetrahedrons-triangles) moved to the right- 
hand side of the picture, for the sake of better visuality. Red oriented line represents 
a -B-type HW configuration, which contains 16 i?-steps, i.e. the maximum possible 
number of them, given by formula (|A.7[) . Remaining 19 steps are of the -B-type, and 
they are depicted as straight lines connecting only two vertices of the traversed G;(5) 
generator, but it is implied that all sites within it are visited. The numbers in the 
figure denote points where HW leaves or enters the lower part of G;+i(5). 



In the remaining two possible situations, when the first and the last Gi(b) are traversed 
by either both 5-steps, or one of them by B-step and the other by £"-step, in a similar 
way it can be shown that inequality (1A.7I) holds as well. 

As one can see in figure lA~3l an B-type HW conformation within Gi+i(b) can use 
all three junctions of any corner Gi(b), which happens when such Gi(b) is traversed by 
B-step. This also means that all Nj junctions within Gi + i(b) can be utilized. Suppose 
that all four corner Gi(b) generators are traversed by B-steps - altogether 12 junctions 
are used by these four E-steps. If the whole HW conformation contains (3 B-steps, the 
remaining (/3 — 4) interior B-steps take 4(/3 — 4) junctions, whereas to all a B-steps 
2a junctions correspond. Every used junction connects two steps, therefore in this case 
inequality 



In a similar manner one can show that in the remaining four possibilities: one, two, three 
or four .B-steps through the corner Gi(b) generators, inequality flA.9j) is valid again. 

Note that formulae (IA.5I) . (IA.7I) and flA.9[) are not correct in the 6 = 2 case. In 
this case all four constitutive Gi(b = 2) generators of the Gi+\(b = 2) generator are at 
its corners, and each of them is connected with each of the other three by two-folded 




(A.8) 



holds, which, again with a + (5 — No, leads to 



P<Nj-NG + 2 = p msx = -(b+l){b + 2)-2 = N G -k E . 



(A.9) 




Figure A3. Example of E-type HW conformation on G;+i(5) with 19 -B-steps, 
which is the maximum possible number /3 max of them, given by formula (|A.9[) . The 
two strands of the conformation are depicted as oriented red and blue lines. The red 
strand starts at vertex labeled with number 1, and terminates at vertex 7, whereas the 
numbers 2 to 6 mark the sites at which the strand leaves or enters the second layer of 
G;(5) generators. 



junction. Since two corner generators have a junction in common, the number of usable 
junctions in relations (A. 3) and (A. 6) is Nj — 2 and Nj — 1 instead of Nj — 4 and Nj — 2, 
respectively. If we represent each G/(2) by a point (vertex) then G; + i(2) represents a 
complete graph of four points. Since there can be neither subgraph with three points of 
degree three and one point of degree one, nor subgraph with two points of degree three 
and two of degree one, terms in relation (12.61) with E 3 and E 2 are not possible. 



Appendix B. Asymptotical behavior of the parameter xi = Bi/Ei 

As it was stressed in section [2] asymptotical behavior of the parameter xi = Bi/Ei is of 
the greatest importance for obtaining the scaling form of HWs. Here we want to prove 
that assumption Pk E +k > ^ks+k+i, for 2 < k < k c + 5, leads to the conclusion that 
x\ — > 0, when I — ► oo. Numbers x\ satisfy the exact relation (12.81) . from which follows 
that 

k +6 kc+5 

~Pk E ~ Pk E +\Xi - pkz+kc+eXi c - £ (pk E +k - m kE+k+ i) x l 

k=2 , n 

X i+1 ~Xi = Xi ^ < , 

J2 Pk E +k xf 
k=0 

where we used the relation = k B — 3, obtained in the previous Appendix, the fact that 
all coefficients p k and m k are positive, as well as the proposed inequality. This means 
that sequence of positive numbers x\ is monotonically decreasing, and since xq = 2, it 
follows that this sequence has non-negative limiting value c = lim^oo xi. If c > 0, again 
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from (12. 8p . when I — > oo, follows equation 



fcc+3 
E m kB+k c k 

c = c 3 ^ 

C fcc+6 



E PkB+kC 
k=0 



which, after dividing by c, can be transformed to 

fcc+5 

= Pk E + Pk E +iC + p kE+kc+6 c kc+6 + (Pk E +k - m kE+k+1 )c k . 

k=2 

However, the latter equation cannot be satisfied, since its right-hand side has the form 
of polynomial in c, with all its coefficients being positive numbers. Consequently, one 
concludes that c = 0. 



Appendix C. Recursion relations for the numbers Ai, Ci and Zi 

In this Appendix we first show that polynomials Rij, appearing in recursion relations 
( 13. 2p for numbers A\ and C\ of one- leg HW conformations, have the form (13.31) . with 
given by (13. 8p . It is obvious that dangling end-step of one-leg configuration within 
Gi + i(b) can be either of A- or C-type. The remaining part of the HW traverses the 
remaining N G — 1 generators Gi(b), and consists of a B-steps and (3 .E-steps, so that 
equation 

a + P = N G -l (C.l) 

is satisfied. 

Consider an A-type HW configuration on generator Gi + i(b). One of its end-steps 
is fixed at one of the four corner generators Gi(b) and it can be either B-, E-, or 
C-step, respectively using one, three, or two junctions. Two of the remaining three 
corner generators are traversed by S-steps, each of them using two junctions of the 
corresponding Gi(b), whereas its third junction cannot be traversed. The last corner 
generator is either traversed by a 5-step, or it contains dangling A- or C-step, and 
such steps use two, one, or three junctions, respectively. If a corner generator contains 
A-step, only one of its corresponding three junctions can be traversed. Now, suppose 
that ,, corner" steps are EBBB, whereas the dangling A-step, using one junction, is in 
the interior generator. Such HW traverses (3 + 2a + 4(/5 — 1) + l)/2 = a + 2f3 junctions, 
and, since corner S-steps cannot use altogether three of all Nj junctions, it follows that 

a + 2(3 < Nj -3, (C.2) 

which with (1C.1I) gives 

(3 < Nj - N G - 2 = l -(b + \){b + 2) - 6 = k c + 2 . (C.3) 

In a similar manner one can analyze all the other possibilities and obtain corresponding 
maximum value of (3. From the table ICll in which all of these possibilities, together 



A model of compact polymers on a family of three-dimensional fractal lattices 22 



Table CI. All possible situations for A-type HW configuration on Gi+i(b), depending 
on the types of its dangling end-step, steps traversing corner generators Gi(b) and end- 
step fixed in one of the corners. For each situation the corresponding numbers of passed 
and usable junctions, as well as the maximum possible value of number (3 of E— steps, 
are given. 



dangling end-step ,, corner" ,, corner" number of number of 

and its position passing-steps fixed-step used junctions usable junctions (3, 



A - interior Gi(b) 


BBB 


E 
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Nj 
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Nj 
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A - interiorG;(&) 


BBB 


B 
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2/3 


Nj 


- 5 


Nj 


-N g 


-4 


A - corner Gi{b) 


BB 


E 




a 


+ 


2/3 


Nj 


-4 


Nj 


-N g 


-3 


A - corner Gi (b) 


BB 


B 




a 


+ 


2/3 


Nj 


- 6 


Nj 


-N g 


- 5 


C - interior Gi(b) 


BBB 


E 


1 


+ 


a 


+ 2/3 


Nj 


- 3 


Nj 


-N g 


-3 


C - interior Gi(b) 


BBB 


B 


1 


+ 


a 


+ 2/3 


Nj 


- 5 


Nj 


-N g 


- 5 


C - corner Gi (b) 


BB 


E 


1 


+ 


a 


+ 2(3 


Nj 


- 2 


Nj 


-N g 


- 2 




C - corner Gi (b) 


BB 


B 


1 


+ 


a 


+ 2/3 


Nj 


-4 


Nj 


-N g 


-4 


C - corner Gi (b) 


BBB 


C 


1 


+ 


a 


+ 2(3 


Nj 


-4 


Nj 


-N g 


-4 



with the corresponding numbers of used and usable junctions, as well as /3 max are given, 
one can see that ku = k\ 2 = kc + 2. 

C-type HW configuration on consists of two strands: one of B-type, i.e. 

with two ends fixed at two corners of the Gi + i(b), and the other of A-type, with one 
end fixed at the third corner, whereas its second end is free, i.e. it can be anywhere 
within the Gi + i(b), in some of interior generators Gi(b), as well as in any of the four 
corner-generators Gi(b). The three steps containing fixed ends can be of either B or 
i?-type, whereas the corner-generator, which does not contain the fixed-end, is either 
traversed by a I?-step or it contains the dangling A- or C-step. The other relevant 
details of possible situations, regarding the type and the position of the step containing 
the dangling end, are listed in table IC2I One can see that HW configurations with A 
dangling end-step and maximal number of .E-steps are those with the dangling end-step 
in the interior Gi(b), whereas the three steps with ends fixed at corner-generators are 
of .E-type, and the fourth corner-generator is traversed by a 5-step, and thus conclude 
that /C21 = Nj — Ng + 1 = kc + 5. Configurations with C dangling end and maximal 
number of i?-steps have also all three fixed-steps of E-type, while the dangling end is 
in the fourth corner-generator, therefore k 2 2 = kc + 5. 

In the remaining part of this Appendix we show that polynomials F and S, 
appearing in the recursion relation (13. ip for the overall number Zf +1 of open HWs 
on Gi + i(b), have the form given by ( 13. 5j) . Polynomials Faa, Fac and Fee correspond 
to HWs whose two step-ends are two A-steps, one A- end one C-step, and two C- 
steps, respectively Since the whole HW configuration has iVc steps, its ,, interior" 
part has iV^ — 2 steps which are either of B- or E-type, implying that Fxy are 
homogeneous polynomials in B\ and E\ of power Ng — 2. In table IC3I we list all 
possible arrangements of steps within open HW configuration with two one-leg end- 
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Table C2. Description of possible C-type HW configurations on Gi + i(b), depending 
on the types of dangling end-step, and the steps within corner-generators Gi(b). For 
each situation the corresponding numbers of passed and usable junctions, as well as 
the maximum possible value of number (3 of -E-steps, are given. 



dangling end-step 


,, corner" 


,, corner" 


number of 


number of 








and its position 


fixed-steps 


passing-step 


used junctions 


usable junctions 


a 


max 




A - interior Gi{b) 


BBB 


B 


a + 2(3 - 1 


Nj - 7 


Nj - 


N G - 


5 


A - interior Gi{b) 


BBE 


B 


a + 2/3-1 


Nj - 5 


Nj - 


N G - 


3 


A - interior Gi(b) 


BEE 


B 


a + 2/3- 1 


Nj - 3 


Nj - 


N G - 


1 


A - interior G; (6) 


E E E 


B 


a + 2(3-1 


Nj - 1 


Nt - 


N r + 


1 


A - corner Gi(b) 


BBB 


A 


a + 2(3-1 


Nj-8 


Nj - 


N G - 


6 


A - corner G;(6) 


BBE 


A 


a + 2/3-1 


Nj - 6 


Nj - 


N G - 


4 


A - corner Gi (b) 


BEE 


A 


a + 2/3-1 


Nj -A 


Nj - 


N G - 


2 


A - corner Gi(b) 


EEE 


A 


a + 2(3-1 


Nj - 2 


Nj 


— N G 




C - interior Gi(b) 


BBB 


B 


a + 2(3 


Nj-7 


Nj- 


N G - 


6 


C - interior Gi (b) 


BBE 


B 


a + 2(3 


Nj-5 


Nj- 


N G - 


4 


C - interior Gi (b) 


BEE 


B 


a + 2(3 


Nj-3 


Nj- 


N G - 


2 


C - interior Gi (b) 


EEE 


B 


a + 2(3 


Nj-1 


Nj 


-N G 




C - corner Gi (6) 


BBB 


C 


a + 2(3 


Nj-6 


Nj- 


N G - 


■5 


C - corner Gi (b) 


BBE 


C 


a + 2(3 


Nj-A 


Nj- 


N G - 


3 


C - corner Gi (b) 


BEE 


c 


a + 2(3 


Nj-2 


Nj- 


N G — 


1 


C - corner Gi (b) 


EEE 


c 


a + 2(3 


Nj 


Nj- 


N G + 


1 


C - corner Gi(b) 


BBC 


B 


a + 2(3 


Nj-6 


Nj- 


N G - 


5 


C - corner Gi(b) 


BEG 


B 


a + 2(3 


Nj-A 


Nj- 


N G - 


3 


C - corner Gi (b) 


EEC 


B 


a + 2(3 


Nj-2 


Nj- 


N G - 


1 



steps, together with the corresponding numbers of used and usable junctions. Using 
these data one can recognize arrangements with maximum possible number of i?-steps, 
and conclude that this number is equal to Nj — N G — 3 = k c + 1. 

If an open HW has both its ends in the same Gi(b), then this corresponds to D- 
or if-step within it, whereas remaining part of that HW consists of a -B-steps and j3 
E'-steps such that a + (3 = No — 1. This means that Sd and Sh are homogeneous 
polynomials in Bi and E\ of power N G — 1. An if -step can exist only in interior Gi(b), 
in which case all steps through corner-generators are of B- type, and four junctions 
can never be utilized. Since such HW takes 2 + a + 2(3 junctions, it follows that 
(3 < Nj — N G — 5 = k G — 1. Finally, -D-step can exist in any Gi(b), either corner 
or interior one, but in both cases the number of used junctions is 1 + a + 2(3, and four 
junctions corresponding to corner-generators cannot be taken by such HW. Thus, one 
obtains that maximum number of .E-steps is /5 max = Nj — N G — 4 = k G . 

Appendix D. Two-leg HW configurations 

Here we give some details of derivation of the conclusion that numbers of two-leg HW 
configurations are not needed for obtaining the scaling form of the overall number of 
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Table C3. Possible steps arrangements of open HW configuration with one-leg ends, 
together with the corresponding numbers of used and usable junctions. The last column 
represents the maximum number of -E-steps for each arrangement, obtained by the fact 
that number of used junctions cannot be larger than the number of usable junctions. 

steps through steps through number of used number of usable /3max 

corner-generators interior generators junctions junctions 



AABB 


(a- 2) B, p E 


1 - 
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-2/3 
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- 6 
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BBBB 


(a - 4) B, P E, 1 A, 1 C 
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Nj 


-4 


Nj 
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-2/3 


Nj 


- 2 


Nj 


-N G 


- 3 


CBBB 


(a -3) B, P E, 1 C 
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ha- 
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Nj 


- 3 


Nj 
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- 4 


BBBB 


(a- 4) B, P E,2 C 


3- 


ha- 


-2/3 


Nj 


-4 


Nj 


-N G 


- 5 



open HWs. 

First we notice that recursion relation for the two-leg HW configuration of the 
type D, given in (13.41) . follows from the fact that dangling ends of the two strands can 
reside either (1) in the same generator Gi(b), thus forming D— or if -step within it, or 
(2) in two different generators Gi(b), thus forming two one-leg steps, either of A- or 
C-type. It is not difficult to see that case (1) can be obtained by cutting a step of 
some .B-type configuration. Therefore, each B-type configuration with k 5-steps and 
(Nq — k) -E-steps gives rise to k D-type configurations with (k — 1) 5-steps, (iV^ — k) 
i?-steps and one D-step, which are obtained by cutting a B-step, whereas by cutting 
its .E-steps (which can be done in two ways for each E— step) one can obtain 2(Ng — k) 
different D-type configurations with k -B-steps, (Nq — k — 1) E-steps and one if -step. 
In a similar way, by cutting an E-type HW configuration, one can obtain an if -type 
configuration. From these observations straightforwardly follows that 

N G N G 

d D (B, E) = 2 km k B k - l E N e- k , d H (B, E) = 2 £ (N G - k)m k B k E N °~ k ~ l , 

k=ks k=ks 
1 Nq N g 

h D (B,E) = -J2 kpf.B k ~ 1 E Na ~ k , h H (B,E) = £ (N G - k)p k B h E N ^ k ' x . (D.l) 

k=kg k=kg 

Applying a reasoning similar to one used in the previous Appendix, one can obtain that 
polynomials dxy and hxY, corresponding to two- leg configurations with dangling ends 
in different generators, have the following forms 

kc+i k X Y 

d XY (B, E) = J2 s^ Y B N o-^E k , h XY (B, E)=Y: ol Y B N o~^E k , (D.2) 

fe=0 k=0 
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with s^ Y and of Y being constant positive integers, XY = AA, AC, CC, and 

kAA = k c + 7 = k A c + 1 = k C c + 2 . 

Next, recursion relations f 1 3 . 1 j) are obtained by dividing (13. 4p with the recursion 
relation for Ei, given in (12.51) . Therefore, coefficients pij are defined as 

d D (Bi, Ei) Ei d H (Bi,Ei)Ei 
Pl1 ^ = F (R F\ ' Pl2 ^ Xl) = F (R F\ ' 

i \ _ ^d(Bi, E t )Ei _ h H (B h Ei) Ei 

P2l{Xl) - — , , , P22{Xl) - — ( s , 

and tji as 

, / \ d X Y(Bi,Ei)Ef h X Y(Bi,Ei)Ef 

hiixi) = — — — , t 2i {xi) - 



Ei+i(Bi,Ei) Ei + i(Bi,Ei) 
where i — 1, 2, 3 corresponds to XY = AA, AC, CC, respectively. Then, from (1D.1I) and 
(fU2l) relations (l3~TH directly follow. 

For I > 1, X| tends to 0, so that from (13.111) one obtains 

Pn(xi) « k B ^x 2 , Pl2 (xi) « 2(A; C + 3)^rf , p 21 ( Xl ) « ^1 , 

xy i 



Pk E Pk E xf ' 

Pk E xf p kE Xi 

Finally, using these approximate forms while multiplying (13. 1Q|) with (|2.9|) . relation 
(I3.12p is obtained, with coefficients a»j equal to 

3 / \ 3 



a n = k B ( — 1 , ai2 = 2(fc c + 3) 



^21 = — , a 22 = [k c + 6) 



and 



>13 



2 \Pk E J " \ Pk L 



* CG fry, \ 2 
b k c +4: / m k B \ 



Pk E \ Pk E j 

Appendix E. Maximally isolated sites within maximally entangled HWs 

In this Appendix we want to derive the formulae (14. ip and (14.4ft for the numbers Nf 1 
and N® 1 of maximally isolated sites within the maximally entangled closed (MEC) 
and maximally entangled open (MEO) Hamiltonian walks, respectively. Since MEC 
HWs on the generator Gi(b) of order / are closed HWs with the maximal number of in- 
steps, and maximally isolated sites within MEC HW are "interior" vertices of 5-steps 
made through unit tetrahedrons (i.e. not the vertices at which HW enters or exits the 
tetrahedron), number Nf 1 is equal to 

Nf I = 2N*, (El) 
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where N B is the number of unit tetrahedrons traversed by B-steps. In order to calculate 
N B , and consequently N l < , we observe maximally entangled B-type (MEB) HWs, i.e. 
B-tjpe HWs with the maximal number of -E-steps, and, similarly, maximally entangled 
-E-type (MEE) HWs. For MEB HW passing through Gi + i(b) generator, we introduce 
label iVj^f for the numbers of unit tetrahedrons traversed by B-step. In a similar way, 
let N^_ B be the number of unit tetrahedrons of traversed by S-step, for MEE 

HW. Then, from the recursion relation for the numbers Bi and Ei (I2.5p . follows the 
matrix recursion relation for the numbers N BB and N EB : 



N B J \ ( k B Na-kn\ ( N BB 



(E.2) 



V Ng» J \k E N G -k E ){ N* 

Since ks and k E are given by ( 12.71) . and for the unit tetrahedron one has N BB = 1 and 
N EB = 0, solving the obtained recursion relation one obtains 

N BB = ^±^3' + NG - kc ~ 6 N l G , (E.3) 

N*» = - NG - kc ~ 6 3 l + NG - kc ~ 6 N G . (E.4) 

' N G -3 N G -3 G K J 

Taking into account that MEC HW on G l+1 {b) consists of k c MEE HWs on Gi(b) 
generators and (N G — k G ) MEB HWs on Gi(b), as implied by f)2.3p . it follows 

N B +1 = {No - k G )N BB + kcNf B , 

and consequently 

N° = ^±^3 l + NG ' k °' 6 N G . (E.5) 

1 N G -3 N G -3 G K ' 

Since Nf 1 = 2N t B , N t = AN l G , and a = as given by (12151) . formula (Oil 

straightforwardly follows. In a similar manner, from ( 13.11) . ( 13.21) . ( 13.31) . ( 13. 51) . ( 13. 8ft . 
and (1K4D for MEO one obtains (IPj) . 
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